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The  paper  deals  with  the  dissipation  mechanism  of  micro  spots  on 
low  pressure  arc  cathodes. 

Based  on  a  new  MFD  stability  theerv  and  on  an  experimentally  fount 
empirical  relation  between  spot  current  and  crater  radius  on  puri 
copper  and  molybdenum  cathodes,  the  time  behavior  of  pressure  and 
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temperature  within  the  spot  and  the  appertaining  voltage  drop 
have  been  calculated.  Emphasis  has  been  placed  on  finding  an 
analytical  solution  for  the  time  dependent  heat  flux  into  the 
hemispherical  crater  surface  as  a  function  of  the  expansion 
of  the  spot. 


This  report  deals  with  the  time  dependent  energy 
dissipation  mechanisms  which  should  occur  in  micro 
spots  on  cold  metal  cathode  surfaces  of  low  pressure 
arcs  . 

The  Appendices  A  and  B  contain  especially  the 
mathematical  way  of  how  the  analytical  solution  of 
the  time  dependent  heat  transfer  equation  and  the 
formula  for  the  time  dependent  heat  flux  for  a 
growing  spherical  spot  has  been  found. 

This  research  has  been  sponsored  in  part  by  the  Air 
Force  Office  of  Scientific  Research/US  Air  Force 


The  paper  deals  with  the  dissipation  mechanism  of 
microspots  on  low  pressure  arc  cathodes. 

Based  on  a  new  MPD  stability  theory  and  on  an  experi¬ 
mentally  found  empirical  relation  between  spot  cur¬ 
rent  and  crater  radius  on  pure  copper  and  molybdenum 
cathodes,  the  time  behavior  of  pressure  and  tempera¬ 
ture  within  the  spot  and  the  appertaining  voltage 
drop  have  been  calculated.  Emphasis  has  been  placed 
on  finding  an  analytical  solution  for  the  time  depen¬ 
dent  heat  flux  into  the  hemispherical  crater  surface 


SURVEY  OF  PICTURES  AND  DIAGRAMS 


Fig.  1:  Hemispherical  crater  model 

2:  Zone  I:  Current  carrying  plasma  channel 
3:  Zone  II:  Hemispherical  crater 
4:  Vapor  pressure  curves 

5:  Zone  II  and  III:  crater  and  molten  layer 
6:  Zone  III  and  IV:  Phase  area  liquid-solid 
7:  Zone  IV:  Solid  metal  cathode 
8:  Temperature  distribution 
9:  Function  y 

10:  Spot  behavior  (cold  copper  cathode) 

11:  Spot  behavior  (cold  molybdenum  cathode) 

A  1:  Temperature  distribution  around  a  ball  with 
the  surface  temperature  Ts 
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Surface  of  the  melting  within  the  crater 
Phase  boundary  area  solid/fluid  within  the 
crater 

Croce  cection  area  of  the  channel  within  the 

plane  of  the  cathode  surface 

Magnetic  induction 

Evaporation  constant 

Vapor  pressure  constants 

Specific  heat  of  the  cathode  material 
Intensity  of  electric  field 
Unit  vector 
Unit  vector 

Factor  depending  on  current  density 

distribution 

Function 

Enthalpy 

Current 

Current  density 
Molar  mass 
Massflow 


mg  Mass  of  an  atom 

4 — * 

P  Pressure  tensor 


Pc  =  Pn 


Vapor  pressure  above  the  melting 


Pressure  within  the  area  Az 

Ambient  pressure 

Heat  of  fusion 

Heat  of  vaporization 

Heat  flow 

Heat  flow  into  the  area  Ac 
Heat  flow  out  of  the  phase  area  As 
Heat  flow  out  of  the  melt  into  the  phase 
area  A  s 

Universal  gas  constant 

Crater  radius  of  the  phase  area  As 

Radius  of  a  sphere 

Velocity  of  the  phase  area  As 

Crater  radius  of  the  phase  area  Ac 

Radiation 

Radiation  into  the  melt 
Temperature 

Evaporation  temperature  =  temperature  of  the 
gas  within  the  crater  straight  above  the  melt 
Temperature  at  the  surface  of  the  sphere  with 

rso 

Melting  temperature 

Cathode  temperature  at  infinity 

Substituted  temperature 

Time 

Voltage  drop  at  the  cathode 


Voltage  drop  in  zone  II 
Voltage  drop  in  zone  III 
Volume 
Velocity 

Exhaust  velocity  of  cathode  material  on  the 
crater  surface 

Velocity  component  in  z-direction  within  the 
area  Az 

Mathematical  variables 
Funct ion 

Matter  constant  of  cathode  material 
Function 

Stability  functions 

Heat  conductivity  of  cathode  material 
Magnetic  permeability  of  free  space 
Thermo  diffusion  constant  of  cathode  material 
Density 

Gas  density  above  the  area  Ac 
Density  of  cathode  material 

Electrical  conductivity  of  cathode  material 
Time 

Potential  function 
Error  f  unct ion 


Work  function 


According  to  todays  knowledge  the  attachment  region 
of  a  low  pressure  arc  or  a  so-called  metal  vapor  arc 
consists  of  one  or  more  bright  "spots"  of  very  high 
current-density  (  10^2  a  /  m  ^  and  even  higher).  A  vapor 
jet  (containing  even  droplets)  and  a  current-carrying 
p 1  a s mac h a n n e 1  emanate  from  each  of  these  microspots. 
The  discharge  channel  which  electrically  connects  the 
cathode  with  the  anode  or  with  the  main  inter- 
electrode  plasma  region  is  therefore  exposed  to  a 
very  strong  axial  flow. 

These  microspots  are  highly  nonstationary,  they 
propagate  and  disappear  and,  depending  on  the 
electrode  material  (oxidized  surface  layer,  surface 
chemistry),  surface  roughness,  overall  temperature, 
and  type  and  pressure  of  gas,  these  microspots  can 
cluster  together  to  macrospots  of  different  sizes 
witn  different  erosion  rates. 

This  report  presents  an  analytical  investigation  of 
the  instationary  energy  dissipation  mechanisms  within 
an  individual  microspot  on  a  pure  copper  and 


HEMISPHERICAL  CRATER  MODEL 


The  following  analysis  considers  a  simplified  model  which, 
however,  according  to  experimental  findings  3 )  ,  24 )  ,  is 
quite  reasonable.  At  first  the  arc  attachment  region  is 
divided  into  five  different  zones  (Fig.  1). 

Zone  I  is  the  current  carrying  plasma  channel  which 
electrically  connects  the  low  pressure  plasma  of  the 
interelectrode  space  with  the  crater  zone  II.  This  latter 
zone  contains  a  complex,  turbulent  mixture  of  evaporated 
material  like  atoms,  ions  and  electrons  under  fairly  high 
pressure.  Adjacent  to  the  high  pressure  zone  II  is  the 
molten  layer  zone  III,  which  is  connected  to  the  solid 
cathode  zone  IV  through  a  hemispherical  phase  transition 
area.  Zone  V  is  the  gas  and  vapor  region  outside  the 
current  carrying  channel  and  is  characterized  by  the  fact 
that  no,  or  only  an  insignificant,  current  density  is 
present  here  . 
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conducting  ambient  zone  V.  As  often  observed  in  experi¬ 
ments, this  jet  can  contain  droplets  which  indicate  that 
the  surface  of  the  melt  within  the  crater  is  boiling  and, 
as  will  be  shown  later, that  due  to  an  explosively 

increasing  pressure  in  the  crater  (zone  II),  molten 

material  is  pressed  and  ejected  out  around  the  crater 

rim.  The  behavior  is  indicated  by  the  protruding  rims  of 

the  craters  left  on  the  surface  of  cold  metal  cathodes 
that  have  been  struck  by  an  arc. 


Zone  I 
Zone  II 
Zone  II 
Zone  IV 
Zone  V 


CURRENT  CARRYING  PLASMA  CHANNEL 
HEMISPHERICAL  CRATER 
MOLTEN  METAL  LAYER 
SOLID  METAL  CATHODE 

NON  CURRENT  CARRYING  GAS  AND  VAPOR  REGION 


Fio.  1: 


Hemispherical  crater  model 


K 


In  our  analysis  we  assume  a  hemispherical  crater,  neglec¬ 
ting  this  crater  rim  and  the  thrown  out  droplets.  However, 
since  these  droplets  have  a  minor  effect  on  the  energy 
balance,  and  since  the  crater  radius  follows  as  a  function 
of  the  spot  current  from  empirical  relations  of  experimen¬ 
tal  observations,  one  indirectly  accounts  for  the  material 
loss  due  to  this  rim  and  droplet  effect  in  zone  III  any- 


ANALYSIS 


Now  we  will  investigate  the  zones  separately  and  match  the 
boundary  conditions  in  between.  Thus  we  obtain  an  integral 
view  of  the  energy  dissipation  and  at  the  same  time  the 
voltage  drop  Uk  between  the  plasma  channel  (zone  I)  and 
the  solid  cathode  (zone  IV).  It  will  be  shown  that  this 
voltage  drop  fits  very  well  with  the  cathode  drop  of  a 
spotty  arc  discharge. 


First  we  look  at  the  current  carrying  plasma  channel, 
zone  I  (Fig.  2).  Former  investigations  in  arc  stability 
5),  6)  lead  to  the  following  criterion  for  a  stable  arc: 


with 
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and  f  -  0 . 8 


which  means  that  the  axial  impulse  transport  within  the 
discharge  channel  J fvfctAi  has  to  be  larger  than 

A* 

7?  / 4+  o.ot-  p-  l 


(3) 


For  a  pure  copper  cathode,  we  know  from  experiments  1  ) 
that  VT  /I  0.12  ( p n't / A )  .  Using  this  we  can  neglect 


the  second  term  within  the  curved  brackets  as  long  as  the 
constant  outer  magnetic  field  is  small  enough  (B0_1  T) 
so  the  criterion  above  is  reduced  to  :  N  0  . 

Utilizing  the  momentum  equation  for  the  crater  zone  II, 
the  vapor  pressure  pc  above  the  molten  layer  can  be  rela¬ 
ted  to  the  axial  mass  flow  (Fig.  3).  Neglecting  the  ambi¬ 
ent  pressure  p  ,  the  criterion  above  can  be  expressed  as 


cc  *  11 

Pc  in  Ai 


which  means  that  in  the  case  of  our  hemispherical  crater 
model  the  vapor  pressure  above  the  melt  has  to  be  about 


3  times  the  average  magnetic  pressure 


If  this  condition  does  not  hold,  a  small  disturbance, 

which  means  a  small  curvature  of  the  channel  axis,  will 
bend  the  channel  more  and  more  until  it  touches  the 

electrode  surface  at  another  point  and  a  new  spot  is 

created.  Fig.  4  shows  the  vapor  pressure  curves  for  copper 
and  molybdenum  with  another  ordinate,  the  product  of 

current  and  current  density,  Ij,  with  dimension  A  2  /  c  m  2  . 
For  Ij  above  these  curves  the  current  carrying  plasma 
channel  is  unstable;  correspondingly  if  Ij  is  less  than 
the  value  given  by  these  curves  we  have  a  stable  discharge 
channel  . 

With  the  assumption  that  the  mean  evaporation  rate  from 
the  fluid  surface  Ac  within  the  crater  corresponds  to  the 
vapor  temperature  Tc,  which  is  the  temperature  of  the  gas 
within  the  crater  straight  above  the  melt  and  is 


the  criterion  (4)  becomes 


M  is  the  atomic  or  molecular  mass  of  the  evaporating 
cathode  material. 


The  minimal  evaporating  rates  for  one  spot  at  a  pure 
copper  and  a  pure  molybdenum  cathode  yield  to 


ho 


4  6 

-9  7* 1 

^  3  •  40 


[*] 
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(6a) 


(6b) 


with  I  in  amperes  and  T  in  Kelvin. 

These  results  are  now  used  for  the  analysis  of  the  dissi¬ 
pation  mechanism  together  with  the  justified  assumption 
that  the  voltage  fraction  of  the  spot  is  just  high  enough 
to  obtain  a  stable  spot. 

First  the  energy  equation  for  the  crater  zone  II  is  formu¬ 
lated,  which  means  that  the  power  converted  within  zone  II 
given  by  the  electrical  power  input 
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'ie.  4:  vapor  pressure  curves  for  copper  and  molybdenum 
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according  to  1 0 )  and  extrapolated  for  pressures 


above  1  bar 


6 


must  be  equal  to  the  difference  of  energy  which  leaves  and 
which  enters  the  crater  zone  per  unit  time. 


Within  this  calculation  the  purely  convective  part 


f  S  h 

Ac 


and  the  parts  due  to  heat  transfer  and  radiation  have  to 
be  separated.  The  convective  energy  flow  can  be  expressed 
as  m!h,  where  m  is  the  evaporated  mass  through  the  crater 
area  A  G  .  Because  of  the  continuity  equation,  m  is  equal  to 
the  massflow  through  the  circular  crater  orifice  area  Az. 
-h  is  the  enthalpy  difference  of  the  plasma-vapor  mixture 
between  the  orifice  area  Az  and  the  region  directly  above 
the  melting  zone  III.  Because  the  mean  pressure  and 
temperature  at  both  areas  are  known,  this  enthalpy 
difference  can  be  calculated.  In  this  evaluation  the 
resultant  heat  conduction  and  radiation  loss  through  the 
area  Az  can  be  neglected  with  respect  to  the  convective 
effects.  The  radiation  and  heat  transfers  (which  include 
ion  bombardment  )  ‘■trough  the  area  A  G  ,  however,  contribute 
importantly  to  the  heating  up  of  the  melting  zone  and  to 
the  evaporation  of  the  cathode  material. 


The  energy  balance  of  zone  III  can  be  written  in  a  similar 
way.  The  electrical  energy  (I  Ujjj)  within  this  zone 
equals  the  difference  of  outflowing  and  inflowing  energy 
per  second  . 
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Two  other  terms  have  to  be  added  in  the  energy  balance  of 
zone  III:  m  Q  v  because  of  the  transition  of  evaporated 
material  to  zone  II.  and  I  $  e  because  of  the  work  function 
of  the  emerging  electrons  (Fig.  5).  The  energy  of  the 
particles  which  flow  over  the  crater  rim  and  that  of  the 
ejected  droplets  are  neglected.  It  is  a  main  fraction  in 
respect  to  the  mass  but  a  minor  fraction  in  respect  to  the 
energy  in  this  zone. 

The  energy  balance  of  the  spotty  cathode  voltage  drop  is 
the  sum  of  the  voltage  drops  of  the  crater  zone  II  and  the 
liquid  zone  III  (Fig.  5).  This  yields  to 

Ui  ( +  A  b  )  +  —  ^  **  -L  ^ ) 

where  (2Tgqs)  is  the  energy  flow  per  second  through  the 
phase  area  solid-fluid  because  of  heat  conduction. 

The  evaporated  mass  m  cannot  become  larger  than  the  molten 
mass  in  the  same  time  period;  this  leads  to  the  require- 
m  e  n  t  : 

**  4  TrJ  rs  ff  no) 

From  Fig.  6  the  energy  which  causes  the  melting  of  the 
solid  material  can  be  calculated  to 

j>f  Qs  •  2irrf  (11) 

where  3S  is  the  heat  of  fusion  and  qp  the  radial  heat  flow 
directed  from  the  phase  area  liquid-solid  to  the  solid 
zone  IV. 


SOLID  CATHODE 


INTO  THE  ENERGY  EQUATION 


I  U,  =  m  (  Ah  +  0V  +  05  )  t  I  j,'  +  ZT  qf 


U'  2  -y-  ( Ah  +  G„  t  Q,j 
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G>  +  IT  C  — 

J  e  S  j 

Fi 

c.  6 

£  Zone  III  and  IV:  phase  ar 

ea  fluid-solid 

With  Eqs.  (10)  and  (11)  and  the  results  from  Fig.  5  ,  the 
heat  flow  in  the  liquid  zone  can  be  expressed  with  respect 
to  the  heat  flow  in  the  solid  zone.  Therefore  the  cathode 
voltage  drop  follows  as: 

2  f  f  6*  )  * $e*  !12> 


Before  looking  at  the  heat  conduction  q^  within  the  solid 
zone  IV,  the  limitation  of  A  caused  by  the  stability  cri¬ 
terion  (6)  and  the  above  requirement  (10)  will  be  added. 
On  a  cold  solid  metal  cathode  the  evaporation  rate  m  with¬ 
in  a  spot  is  limited  by 


£  /i Trj 


0 


ft 


(13) 


To  calculate  the  heat  flow  q  f ,  the  heat  flow  equation 
within  zone  IV  (Fig.  7)  has  to  be  solved.  Therefore  two 
assumptions  are  made: 

a)  the  heat  transfer  coefficient  '  f  respectively  the  ther¬ 
mal  diffusion  coefficient  X  f  within  the  solid  is  constant 
and 

b)  the  term  which  takes  the  ohmic  heating  within  the 
solid  into  account,  can  be  approximated  step  by  step  by  a 
potential  function  with  f'-l/r. 


Hence  it  is 


in  the  region 


the  potential  function  is 
given  as 


5 

i 


« 


c 


r  -  * 

V  (s,  0  -  y5^  (rt  t) 

n,  r  £■  rt 

f  * 

f  (r>  t)  *  ft  ( r,  t) 

$  • 

*  * 

^  *  f*v 

I 

I 


Setting 
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*  r  IrT1 


f-f 


C***(f. 


there  follows  because 


,  iz  d 
A  rv 


(16) 


which  is  in  fact  no  potential  function.  But  within  a  small 
region  *  r  k  the  term  (14)  can  be  approximated  by 


where  *  O  is  nearly  realized.  This  can  be  done  within 
all  regions  v  =  1,2,3...  •  With  this  step  by  step 
approximation  of  £  f  j  2  by  a  potential  function  'f  the  heat 
flow  equation  can  be  rearranged. 

Replacing  T  by 


(  17  ) 


where  T  *  is  the  new  independent  variable 
the  new  differential  equation  follows  as 


£ 

S>T 


(  1 
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which  is  solved  for  the  given  initial  and  boundary  con¬ 
ditions  in  Appendix  A.  Equation  (17)  leads  to  the  solution 
for  the  temperature  T(r,t)  which  is  shown  in  Fig.  8.  Using 
this  solution  T(r,t)  the  heat  conduction  loss  can  be  cal¬ 
culated  by  forming  suitable  gradients  on  the  solid  side  of 
the  phase  area  fluid-solid  (r=rs;  exactly: 

lr—0 


r  =  1 i m  (rs  +  ir)  . 


In  Appendix  B  the  heat  flow  into  the  solid  zone  IV  is  cal¬ 
culated  taking  care  of  the  fact  that  at  time  t=0  the 


initial  radius  is  r 


s  o 


=  0; 


x  (p-Tj  fr*  I  (  -„)] 


(  19  ) 
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^5*  77*T 
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The  function  /(  '  is  displayed  in  Fig.  9.  Taking 

the  rises  of  current  I(t)  and  spot  radius  rs(t)  from  ex¬ 
periments  2 )  Eqs.  \‘~i )  allows  the  calculation  of  the  time 
dependent  heat  flow  into  the  solid  electrode.  With  this  it 
is  possible  to  determine  the  cathode  drop  voltage 
(without  inductive  part),  the  minimal  vapor  pressure  pc 
and  the  spot  temperature  Tc  as  a  function  of  time  using 
Eqs.  (12), (4)  and  vapor  pressure  diagrams  (Fig.  4). 


RESULTS  AND  DISCUSSION 


It  is  known  from  experimental  investigations  1  )  ,  2  )  ,  1  1  ) 
that  during  the  first  3  to  5  ns  of  a  pulsed  discharge  with 
a  current  rise  of  about  1  to  5  *109  A/s,  a  small  circular 
melting  point  occurs  on  the  pure  metal  electrodes  which 
afterwards  forms  a  crater;  then  the  crater  rim  expands 
with  a  velocity  of  10^  to  5*10^  cm/s.  After  about  5  to  40 
ns  the  final  crater  size  is  reached  which  means  that  the 
current  transportation  is  incured  by  a  new  spot. 

On  pure  metal  cathodes  the  crater  radius  rs  increases  with 
the  current  I  s  •  Experiments  lead  to  the  empirical 
approach  : 

0  *  'i  e 

with  pure  copper:  r0=1,7  m  and  I*  =  83  A 
with  pure  molybdenum:  r 0  =  1  . 2  pm  and  I#  =  33  A  2  )  . 

It  is  now  possible  to  analyze  the  time  behavior  of  the 
dissipation  mechanism  of  a  cathode  spot  and  to  calculate 
the  purely  ohmic  part  of  the  cathode  voltage  as  a  function 
of  time  if  one  takes  the  spot  current  Is  as  a  given  quan¬ 
tity  which  linearly  rises  with  time. 

Figures  10  and  11  show  the  results  of  these  calculations 
for  copper  and  molybdenum.  In  both  cases  the  characte¬ 
ristic  traces  of  the  curves  are  similar  to  each  other  for 
the  vapor  pressure  pc  ,  the  surface  temperature  Tc  and  for 
the  ohmic  part  of  the  voltage  U ^ . 


During  the  initial  phase  the  pressure  rises  explosively 
and  sooner  or  later,  depending  on  electrode  material  and 
current  rise  (dl/dt),  it  reaches  a  maximum.  With  copper 


(dl/dT  =  2.67  1C)9  A/s)  this  maximum  of  about  160  bar  is 
reached  after  about  30  ns;  with  molybdenum  (dl/dt  = 
8  -  1 0  9  A/s)  it  needs  only  about  3-7  ns  to  get  to  the  maxi¬ 
mum  of  about  55  bar.  This  very  significant,  explosive 
pressure  rise  on  the  inner  crater  surface  gives  an  expla¬ 
nation  for  the  fluid  metal  which  is  ejected  and  which 
forms  a  crown  like  rim. 


Initially  the  temperature  rises  steeply,  later  on  it  rises 
more  slowly  up  to  nearly  5000  K  for  copper,  and  nearly 
7000  K  for  molybdenum.  Nevertheless  the  temperatures  do 


not  reach  the  value  which  would  be  needed  for 
thermionic  emission  even  in  the  case  of  molybdenum. 


pure 


To  explain  the  very  high  electron  emission  on  the  cathode 
G.  Ecker  postulated  a  combination  of  thermionic  and  field 
emissions  12).  This  seems  to  be  possible  because  new 
experiments  1 3 )  show  field  amplifying  "Taylor  cones"  at 
formerly  plane  fluid  surfaces.  With  this  explanation  the 
emission  of  the  charged  particles  is  not  uniformly 
distributed  on  the  inner  fluid  crater  surface  but  will 
occur  at  one  or  more  "Taylor  cones"  which  emerge  and 
vanish  at  the  fluid  surface.  Therefore  only  the  time 
average  of  the  emission  is  uniform  over  the  inner  crater 
surface  . 


A  complex,  turbulent  high  pressure  plasma  exists  in  the 
crater-zone  where  the  Debye-length  reaches  quite  possibly 
only  10~®  to  10“7m  and  therefore  the  electrical  field  at 
the  tips  of  the  "Taylor  cones"  could  be  some  1 0  9  V/m. 
With  that,  the  relatively  high  emissions  of  10®  A/cm2 
(and  more)  can  be  explained.  A  further  possibility  would 


be  the  evaporation  of  negative  ions  out  of  the  tip  of  such 
Taylor  cones,  similar  to  those  observed  with  positive  ions 
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on  anode  tips  1  ^  ^ .  The  life  span  of  such  negative  ions 
within  the  crater  plasma  is  expected  to  be  short;  they  de¬ 
compose  rapidly  into  free  electrons  and  atoms;  through 
further  impact  the  atoms  are  ionized  into  electrons  and 
positive  ions.  These  considerations  are  only  speculative, 
but  they  could  yield  a  plausible  explanation  of  the 
experimentally  found  high  current  densities  of  about  10® 
A/cm2  and  greater. 

The  purely  ohmic  part  of  the  cathode  voltage  above  the 
spot,  first  drops  ste  ply  from  relatively  high  values, 
passes  through  a  minimum  and  then  increases  slightly.  It 
is  of  interest,  that  the  high  voltages  at  the  start  are 
caused  mainly  by  the  high  heat-conductivity  losses  into 
the  electrode  material  at  the  beginning.  In  the  first  mo¬ 
ment  after  the  ignition  of  the  spot  discharge,  the  tempe¬ 
rature  decay  at  the  phase  boundary  layer  solid-liquid  is 
very  steep  (see  Fig.  8)  and  therefore  the  heat  flux  into 
the  solid  cathode  is  very  high.  The  actual  ignition  is  not 
described  by  this  analysis.  The  ignition  voltage  however 
has  to  be  surely  greater  than  or  at  least  equal  to  the 
purely  ohmic  cathode  voltage  U^,  calculated  here,  to  init¬ 
iate  a  spot  discharge  which  leads  to  a  craterlike  cavity. 

Also,  the  time  when  the  channel  becomes  unstable  is  not 
included  in  this  evaluation  -  the  assumption  was  that  the 
stability  criterion  (1)  or  (4)  is  just  satisfied  at  any 
time.  This  means  that  the  product  of  current  and  current 
density  Ij  would  reach  its  maximum  at  the  same  time  as  the 
vapor  pressure  pc  .  However,  it  is  improbable  that  with  a 
further  increase  of  temperature,  current  and  voltage,  the 
emission  ^r  the  average  current  density  j  at  the  liquid 
surface  of  the  crater  would  decrease.  An  increase  of 
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with  the  conditions: 
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Solution : 

-  .  .r 

-  The  term  :  j-  with  j  s  r ~ 'j  is  represented  within 
hemispherical  shell  by  the  potential  function 


fij2-  f 


-  Substitution: 


T(r,i)  =  T*M+  ffijW 


00  n  o  w 


*  F‘>r[rlSlAjV]  ***  *  ° 


1 1°  h  e  a  *■  conduction  equation  becomes 
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within  the  tin0  interval  the  conditions  are 
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The  solution  of  this  problem  is  found  with  Landau-Lif- 
schitzIT)  as 


_ 4 


-* 


rj.sdz*)  I  <  e^jT  (i5) 


tr-GF  I  (t-Tj** 
T"  -oo 


where 
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T 
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With  this  approach  the  temperature  field  is  determined  if 
the  temperature  T*(rso,~)  on  the  surface  of  the  sphere  is 
known  at  any  time.  The  location  of  the  temperature  front, 
T  s ,  at  any  time  is  known  by  experiment.  Therefore  we  try 
to  reduce  the  time  dependent  temperature  behavior  on  the 
surface  of  the  sphere  to  the  behavior  of  the  temperature 


front  T  s  . 


Choosing 


r 

ATSt>(l)  -  f  c/t 

a _ 


the  solution  with 
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( A  6  ) 


( A  7  ) 


results  in 
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where  $(x)  e  J  B  t*y  is 


the  error  function. 


With  the  substitution 


r'rJo 


r-n. 


(  A  1  0  ) 


^  yr  (A9) 

and  using  ( A7  ) 

r'r h)  -*('!--) 

the  following  solution  is  obtained: 
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(All) 


As  can  be  shown  analytically  with  the  linear  approach 
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(All)  are  the 

solutions  of 

the  problem 

same 

way  , 

,  at  any  time 

r  +&  t 

small),  the 
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) 
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(  A  1  3  ) 


This  yields  also  the  equation  ( A  8  )  or  (All). 
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Substituting  f(T)  in  Eq .  A  (  1  1  )  with  respect  to  A  (  1 3  )  we 
obtain 

to 

T(r,i)-TK  *  y(Ts- *. ] 
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_  t  *  (AIM) 

4  +  I f1 j'Mdi 


or  by  using 


and  calculating  •  (l  )  with  (A10) 


where 


r(r,i)-7„  .  $  (Tt  ] 

*IH)U  f  '<*)] 


®°  t  i 

<f(*)  *  jp  j(/f~'p‘)e  cty  *  [j~  &(*) ](4+2xl) 


( A  1  5  ) 


l  -* 

rn  X  C 


(  A  1  6  ) 


"  z/*T 


and 
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[  ?  /  ) 

k0  and  J  z»  have  to  be  chosen  at  any  time  with  respect 
to  the  condition  that  the  location  of  the  temperature 
front  at  time  t  in  the  surrounding  medium  is  the  experi¬ 
mentally  given  one.  So  we  obtain  from 
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W  / 


(*) 
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and  with  (A15) 
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with 
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Following  from  ( A  17 )  ,  because  of  t  -*  0,  rs  =  rso, 
xs  -0, 


and 


that  is, 
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(A  19 


In  the  same  way  we  obtain  with  ( A  1 7  )  and  ( A  1 9  ) 
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Using  this  together  with  ( A  1 9  )  ,  E  q .  (A15)  yields  the 

temperature  d i s t r i b u t i on ,  at  any  time,  and  at  any  place 
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APPENDIX  B 


Calculation  of  the  heat  flow  2  ”  r  g  q  f 


The  heat  flow  at  time  t  through  the  surface  of  a  half 
sphere  with  radius  r(>rs)  into  the  colder  solid  body 
leads  to 


£  T(r,t) 


( B  1  ) 


The  heat  flow  is  obtained  as  a  function  of  t  with  the 
passage  to  the  limit: 
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—  T(r.t)  can  be  calculated  with  Eq.  (A21): 
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In  addition  is 
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